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Abstract 

The stability of shear flows of electrically conducting fluids, with respect to finite amplitude 
three-dimensional localized disturbances is considered. The time evolution of the fluid impulse 
integral, characterizing such disturbances, for the case of low magnetic Reynolds number is obtained 
by integrating analytically the vorticity equation. Analysis of the resulted equation reveals a new 
instability criterion. 

1 Introduction 

The use of magnetohydrodynamics (MHD) for transition and turbulent control is quite attrac- 
tive from the point of view of applications, in particular when new technologies may allow direct 
turbulent control in sea water. Most of the work in this field has been concerned with flows of 
electrically conducting fluids subjected to transverse magnetic fields. This configuration is used in 
MHD generators, accelerators and pumps. Recently, Nosenchuck and Brown (1993), demonstrated 
experimentally that the application of wall-normal Lorentz force prohibit lift-up and bursting of 
near wall fluid, which are characteristics of the end-stage of transition and near wall turbulence. 
However, in such a configuration, the magnetic effects are due mainly to the coupling between the 
mean velocity profile and the magnetic field, rather than the damping of turbulence. 

On the other hand, when the mean flow and the magnetic field are aligned, the direct effect of 
the magnetic field is on the disturbed velocity field. Fraim and Heiser (1968) studied experimentally 
the effect of a strong longitudinal magnetic field on the flow of mercury in a circular tube. They 
found that the magnetic field can significantly increase the critical Reynolds number (Re cr ) for 
transition. More examples showing similar results are summarized in the book of Branover (1978), 
who concluded that linear theories of stability of MHD flows yield values of Re cr which are much 



higher than the measured values. This suggests that a proper explanation of the above mentioned 
experimental results must take nonlinear effects into account. 

The evolution of a finite amplitude three-dimensional localized vortex disturbance embedded 
in an external incompressible shear flow was considered by Levinski & Cohen (f995, hereinafter 
reffered to as LC). Using the fluid impulse as an integral characteristic of such a disturbance, 
they found that parallel shear flows are always unstable with respect to localized disturbances, the 
typical dimension of which 5, is much smaller than a dimensional length scale A, corresponding to 
an 0(1) change of the external velocity. Moreover, their analysis predicts that the growing vortex 
disturbance is inclined at 45° to the external flow direction, in a plane normal to the transverse 
axis. It was also shown that although viscosity plays a crucial role in the generation of the initial 
localized disturbance and in determining the mean flow field, it plays no role in the time evolution 
of its fluid impulse integral. In other words, once the mean field is established, the subsequent 
evolution of the disturbance fluid impulse integral is largely an inviscid one. These predictions 
agree with previous experimental observations concerning the growth of hairpin vortices in laminar 
and turbulent boundary layers, see e.g. Head & Bandyopadhyay (1981), Acarlar & Smith (1987) 
and Hagen & Kurosaka (1993). The application of this approach to Taylor-Couette flow revealed 
a new instability criterion, which was recently verified experimentally by Cohen et al. (1996). 

The purpose of the work reported here is to examine the effect of an externally imposed magnetic 
field on the onset of such disturbances. The analysis is restricted to incompressible shear flows 
characterized by a low magnetic Reynolds number, Re m = [icrAU <C 1, where fi is the magnetic 
permeability, a is the electrical conductivity and U is a characteristic velocity scale of the external 
flow. The magnetic Reynolds number can be considered as the ratio between the diffusion time 
scale of the magnetic field //oA 2 , and the hydrodynamic time scale A/U . For small values of Re m , 
the magnetic field induced by a disturbance diffuses rapidly. This leads to the dissipation of the 
disturbance kinetic energy and consequently may stabilize the flow. 

2 Analysis 

For Re m <C 1 and stationary external magnetic field, it was shown by Braginskii (1960) that 
the electromagnetic force per unit volume is f = a [(Ur x B) x B — V<Fr x B], where B is the 
magnetic induction of the external field and Ut is the total velocity vector as defined below. The 
scalar potential $t is determined from the condition that the charge is neutralized, which for a 
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uniform liquid yields 

V 2 <J> T = B fi T , (1) 

where fl-r = V x Ut is the total vorticity vector. For this case, the three-dimensional vorticity 
equation for an incompressible flow is given by 

_L + (u T • V)ft T - (O t • V)U T - -(B • V) (U T x B - V$ T ) = vAfl T , (2) 
at p 

where p and v are the density and kinematic viscosity of the fluid, respectively. 

We consider the flow field as being the sum of two contributions: the external shear field in which 
fi = V x U, and a finite amplitude disturbed field in which u = Vxu, Thus, the total velocity and 
vorticity vectors can be written as Ut = U + u and Ox = fi + cj, where the undisturbed external 
flow field is assumed to be a known solution of ([!]) and (||) when u = 0. The initial disturbed 
vorticity <*j(x, to) = wq(x) is assumed to be confined to the small region of order 5 -C A as well 
as 5 <C Ai, where Ai is a typical length scale corresponding to an 0(1) change of the external 
magnetic field. 

Owing to the smallness of the disturbed region, the external velocity, vorticity and magnetic 
fields are approximated by Taylor series expansions. Following LC, we use a Galilean frame, mov- 
ing with the disturbance, i.e., U(0) = 0, and consider the initially embedded vorticity region as 
surrounded by an infinite field having a constant velocity shear and a constant magnetic induction. 
Consequently, 

^( X ) = E^^'> fii(x)=«*(0), Bi(x)=Bi(0) where i = 1,2,3. (3) 

3=1 J 

When the equations for the undisturbed external flow and magnetic field are subtracted from (|T|) 
and (H|) respectively, we obtain 

3u> (j 

— + (U • V)u> - (o;-V)U - (fl • V)u + (u • VW - (cj-V)u - -(B • V) (u x B - Vcfi) = 0, (4) 
at p 

and 

V 2 <P = B-lu. (5) 
In (El) the viscous term is omitted since, as was shown by LC, the corresponding contribution of 
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the viscous term to the dynamics of the fluid impulse integral vanishes in view of the asymptotic 
behavior of the disturbance vorticity far from the origin. 

We shall follow the evolution in time of the fluid impulse integral of the disturbance, defined as 

dP 1 f c^(x,i) 

— = - / x x — -dV , (6) 

dt 2 J at w 

where x is the position vector, dV is a volume element and the time derivative of u>(x, t) is 
determined from (||). 

Since the time evolution of the fluid impulse is an integral over the whole volume, we must 
first verify that most of the contribution to this integral comes from the localized disturbed region. 
Indeed, all of the vorticity contributing to this integral, except for the part generated via the fourth 
and the seventh (electromagnetic) terms in (Q), is confined to the disturbed region. 

In order to estimate the contributions of the fourth and the seventh terms in (||), we examine 
their far-field behavior. The expression for the scalar potential (j) is obtained from the general 
solution of (H) 

1 ^ B - W(X1 W 1 = B.M, (7) 



4-7T J |x — Xi| 3 

where the asymptotic series of M, expressed in terms of the fluid impulse, is given (Batchelor, 
1967) by 



M = ipx^l0d). (8) 
47r xr xr 



Similarly, the far-field velocity induced by the localized vortex disturbance is 



u(x) ~ ~h 



P 3(P • x)x 

I |3 l |5 

X X 



WAV (9) 



Substitution of @ and (|9|) into (Q) shows that the far-field vorticity diminishes in magnitude as 
|x|~ 4 . Consequently, the integral (^) is not absolutely convergent and depends on the way in which 
the volume of integration is allowed to tend to infinity. 

In order to overcome this difficulty we follow the procedure proposed in LC. Accordingly, we 
subdivide the velocity and vorticity fields into two parts, u = oj 1 + u> n and u = u 1 + u 11 , so that 
u 1 ' = V x u 1, . Therefore, for each part we require that 

V • u 1 = V • us 11 = 0. (10) 
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The first part, indicated by the superscript I, is associated with the concentrated vorticity confined 
within and in the vicinity of the initially disturbed region, whereas the second part, indicated by 
the superscript II, is associated with the far-field vorticity generated by the problematic terms (the 
fourth and seventh) in (B). Accordingly, we set the initial distribution of the vorticity fields as: 



u/(x, t = to) = o; (x) and tj /7 (x, t = t ) = 0, 



(11) 



and follow the evolution of w 7 (x, i). In addition, we subdivide the whole space into two regions, 
inside and outside a spherical domain of radius R > 8, enclosing the disturbance. In the outer 
region the corresponding system of the vorticity equations is given by 



du 1 
~~dt 



+ (U • V)u/ - (oj j -V)U + (u • V)« - (wV)u - (fl • V)(u 7 - u ) 



a 
P 

duj 11 



dt 



(B • V) [(u 1 - u ) x B - V (B • (M 7 - M ) 
+ (U • V)u H - {u H -V)\J - (O • V)(u 7/ + u ) 



a 
-( 

P 



(B • V) [(u 11 + u ) x B - V (B • (M 11 + M 



where uq and Mq are the leading terms in the asymptotic series of (H) and (ft), 



u 



1 

47T 



P 3(p • x)x 

x| 3 Ixl 5 



Mn 



l x 

— p x ris-- 

47r xr 



and the fluid impulse p corresponds only to u 1 , i.e., 



p = - J xx cj 7 (x)rfF. 



(12) 



(13) 



(14) 



(15) 



Since uo and Mo cancel the leading terms of the far-field vorticity, generated via the problematic 
terms in (|i~2|), the asymptotic behavior of u 1 is given by 



u/(x,t) 



O 



for |x| » 5, 



(16) 



and therefore the fluid impulse integral ([D]) is absolutely convergent. 
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For the inner region we write 



du 1 
~dt 



+ (U • V)^ 7 - (cj / -V)U + (u • V)w - (w-V)u - (O • V)u 7 



-(B • V) (u 7 x B - V(B • M 7 )) + W = 0, 



(17) 
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dt 



+ (U • X7)u n - (u/ 7 -V)U - (n • V)u 77 - -(B • V) f u 77 x B - V(B • M 77 )) - W = 0, (18) 

P v ' 



so that the sum of the two equations in each region yields Eq. (0) for that region, and together 



with the initial conditions given in (11), yields the original problem for the entire space. 

For the outer region, condition (|i"o| ) is always satisfied. For this condition to be satisfied in the 
inner region V 2 ^ must be equal to zero, as can be shown by applying the operator (V-) to ([17]) and 
(|l~8|). Then, ^ is determined by solving the Neumann problem for which the normal derivative of 
at |x| = R is matched with the scaler product of the unit vector normal to the boundary surface 
n, and the terms in (|i~2|), containing uq and Mq, i.e. 



n 



dn 



n 



|x|=B 



a 



(ft ■ V)u + -(B • V) (u x B - V(B • M )) 
P 



(19) 



|x|=i? 



Accordingly, using (|l4|) the expression for ^ is given by 



87T.R 4 



(fl • p)|x| 2 - 3(n • x)(p • x) + — (B • x)(x • (p x B)) 



(20) 



Since the integral ( |15|) is absolutely convergent, it is convenient to use an infinite sphere as the 
volume of integration. Consequently, the time evolution of p is given by 



* = I , im / x x a-'fr'W, 

dt 2 Ri->co J\x\<Ri Ot 



(21) 



Substitution of (12) and (fij) into <M) yields 



dp 1 

— = lim / xx 

dt 2 fli->oo J\x.\<R-i 



(U • V)u; 7 - (^ 7 -V)U - (O • V)u : + (u • V)w - (wV)u 



dV 



H lim / x - 

2p Ki-»oo y| x |<fli 



(B • V) (u 7 x B - V(B • M 7 ; 



dV - - lim / xx [(n • V)u 

2 Rl-^OO J R<\ X \<R % 
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+-(B-V)u xB dV + — lim / x x (B • V)V(B • M Q )dV - - [ xxV$rf7. (22) 

p J 2/9 Ri-»oo Jr<\x\<R-i 2 7|x|<fl 

Each one of the integrals in (p5) is evaluated for a finite value of Ri and its limit as R\ — ► oo is 



then taken. The result of the first integral in (22) was already obtained in LC. As shown in the 
Appendix for such domains of integration, the integral contribution of the last three integrals in 
d22| ) is identically zero. Therefore, the artificial vorticity field has no direct impact on the evolution 
of the fluid impulse associated with the concentrated vorticity oj (x, t). Moreover, as was shown in 
LC, the influence of the vorticity field u n (~K,t) on the evolution of w (x, t) can be neglected. The 
second integral in (|2^) is calculated using a similar procedure to that described in the Appendix 
and in LC. Finally, (p2[) becomes 



* = _> VCp . u) _i (p .v)o-^p + icB-p)a (23) 



3 Application to a representative example 



In the following we consider a simple example in which ( |23|) is applied to a parallel shear flow of elec- 
trically conducting fluids and a new instability criterion for finite amplitude localized disturbances 
is found. For a parallel plane shear flow, the external velocity field is given by U = (U(y),Q, 0) , 
for which a right-handed coordinate system is used with x = (x, y, z), where the vector entries are 
the downstream, cross-flow and spanwise directions, respectively. The direction of the magnetic 
induction vector is chosen to be parallel to the downstream direction, i.e. B = (.6,0,0). As was 
mentioned above, such a flow does not exhibit a direct coupling between the mean flow and the 
magnetic field. Consequently, the direct effect of the magnetic field on the turbulent structure can 



be revealed. In this case, equation ( |23| ) for the fluid impulse vector p = (p x ,Py,Pz), is reduced to 



dp x _1 dU__o^_ dpy 1 dU 2aB 2 dp z 2aB 2 
~dt ~ ~2 Py ~dy~ ~ ~~5p~ Px ' ~dt ~ ~2 Px ~dy~ ~ bp Py ' ~dt ~ l~p~ 



for which the eigenvalues {Aj}f =1 can be found from the characteristic equation 
I 2aB 2 



x2 3(7^2 / a B 2 \ 2 fldUY 

A? + ^— A, + 2 — - - — 



hp \ hp J \2 dy J 



0. (25) 
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Hence, the flow under investigation is stable with respect to three-dimensional localized disturbances 
only if the real part of \ is not positive. Therefore, for stability we require that 



where N is a dimensionless interaction parameter which represents the ratio between the electro- 
magnetic and the inertia forces. 

The term 'stability' (and 'instability') here means that the fluid impulse of a closed localized 
(in all three directions) vortex disturbance will not grow (or grow) in time. The fluid impulse is 
a very suitable characteristic of localized vortex structures such as hairpin vortices in boundary 
layers, since it combines the geometrical dimensions of the structure together with the magnitude 
of its vorticity field. As such, this stability definition cannot describe 'wavy' disturbances or 'quasi' 
two-dimensional structures for which the fluid impulse integral is not defined. It should be noted 
however, that the above stability definition is not equivalent to the conventional criteria of linear 
stability and energy stability. In fact, the growth of the fluid impulse does not necessarily guaranty 
growth in energy of the localized disturbance. For example, viscous diffusion leads to the decay of 
the localized disturbance energy while its fluid impulse remains the same. 

As an example of the new stability criterion, the value of the interaction parameter required 
for stability of Poiseuille flow in a tube is aB 2 D/pU > 14.1, where D is the tube diameter and 
U is the mean velocity. For flow of mercury in a circular tube subjected to a strong longitudinal 
magnetic field, Fraim & Heiser reported an increase of the critical Reynolds number for transition 
from 2250 to 10350 when the magnetic induction B was increased from zero to 1.75weber /m 2 , and 
the corresponding interaction parameter at the upper limit was 9. As can be seen in Fig. 10 of 
their article, this value of the interaction parameter is an order of magnitude larger than the value 
predicted by Stuart's linear theory (Stuart, 1954). A direct comparison of these results with our 
prediction is questionable because of at least two reasons. The first is that the experimental mean 
velocity profiles were not reported by Fraim & Heiser. The second is that the criterion presented 
here is applicable only with respect to localized disturbances whereas transitional flows include 
various types of disturbances. Nevertheless, the predicted value of the interaction parameter found 
in this paper is of the same order of magnitude as that of the experimental ones for large Reynolds 
numbers. 




or N = 



pdU/dy 2^/2 




(26) 
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Appendix 

In this Appendix we show that the last three integrals of (^), defined as J 1 , J 2 and J 3 , respectively, 
are identically zero. In Cartesian tensor notation, the sum of J 2 and J 3 is 



a 



[Z I 



d 

R<|x|<i?i 3 dx k 



1 r 

(B-V)(B-M )dV --e ljk xj — dV, (Al) 

2 J\xl\<R OXk 



where e,ij k is the alternating tensor and the usual summation convention is applied. 
Integration by parts and using Gauss' divergence theorem yields 



J 4 ~\~ J 4 



a 



I n kXj {B ■ V)(B • Mo)dS - I n kXj (B ■ V)(B • M )dS 

J|x| = /?! J|xl=i? 



|x|=R 



o jk i (B-V)(B-M )dV 

IR<\x\<Ri 



', e ijk 



<p n k Xj^dS - 5 jk / ^dV 

J\x\=R J\xl\<R 



(A2) 



where 6{j is the Kronecker delta function. Using the properties of the alternating and symmetrical 
tensors e{j k 5j k = 0, and eij k n k rij = 0, it follows that J 2 = J 3 = 0. 
A similar procedure for J 1 gives 



2 Wlx^il! 



x=R/ 



[dS(Sl ■ M ) - (dS • O)(M + x x u ) 



a r 
+- 
P 



B 2 dS x M - B((dS x Mo) ■ B) - (dS • B)(x x (u x B)) 



(A3) 



When ( |l4| ) is substituted into (A3), the result of the two surface integrals become independent of 
the surface radii and consequently cancel each other, i.e., J 1 = 0. 
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